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Abstract 

The stability of Marangoni roll convection in a liquid-gas system with deformable in- 
terface is studied in the case when there is a nonlinear interaction between two modes of 
Marangoni instability: long-scale surface deformations and short-scale convection. Within 
the framework of a model derived in , it is shown that the nonlinear interaction between 
the two modes substantially changes the width of the band of stable wave numbers of the 
short-scale convection pattern as well as the type of the instability limiting the band. De- 
pending on the parameters of the system, the instability can be either long- or short-wave, 
either monotonic or oscillatory. The stability boundaries strongly differ from the standard 
ones and sometimes exclude the band center. The long-wave limit of the side-band insta- 
bility is studied in detail within the framework of the phase approximation. It is shown 



that the monotonic instabihty is always subcritical, while the long-wave oscillatory insta- 
bility can be supercritical, leading to the formation of either travelling or standing waves 
modulating the pattern wavelength and the surface deformation. 
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1 Introduction. 



An important class of instabilities of one-dimensional patterns are side-band instabilities. The 
most important one among them is the Eckhaus instability |2[ which is a monotonic long-wave 
instability that renders the pattern unstable with respect to infinitesimal compressions and 
dilatations. Its onset can be characterized by the vanishing phase diffusion coefficient |^]. Near 
the threshold, where the pattern can be described by a real (dissipative) Ginzburg-Landau 
(GL) equation, the Eckhaus instability arises generically if the deviation q of the wave-number 
of the pattern from that corresponding to the bifurcation point exceeds a certain value, i.e. if 

\q\ > e/Vs, (1) 

where is the deviation of the control parameter from the threshold. A similar instability 
of wave patterns (e.g. described by the complex GL equation) is usually referred to as the 
Benjamin- Feir instability [Q]. Typically, the Eckhaus instability restricts the range of allowed 
wave numbers of the pattern, i.e. its nonlinear evolution eventually brings the system back to 
the stable range of wave numbers. 

The Eckhaus instability of periodic patterns has been studied extensively, both theoretically 
and experimentally, in various pattern-forming systems, such as Rayleigh-Benard convection 
^ , convection in binary mixtures , Taylor vortex flow ||8|, |9| , electrohydrodynamic convection 
in nematics [^], directional solidification [11|. Even quite close to the threshold, the boundaries 



of the Eckhaus instability of stationary periodic patterns can differ from (|l|) due to various 
factors, such as the system confinement ^ subcritical character of the instability |13], the 
shape of the neutral stability curve 0, |l5|, the effect of nonlinearity |16] and resonances Q, 



the presence of random fluctuations [|17|, chiral symmetry breaking (e.g. due to rotation) |ll8| , 
etc. 

Although the Eckhaus instability is the most common one-dimensional side-band instability, 
others are possible as well. In the present paper we show that spatially periodic patterns arising 
due to Marangoni convection in a liquid layer with a deformable free surface heated from below|^ 
can possess additional side-band instabilities: a short-scale steady side-band instability and a 
long-scale oscillatory side-band instability. An important feature of this system is the presence 
of two primary modes of instability (generating the Marangoni convection itself): a long-scale 
instability (with zero wave-number) connected with the deformation of the free surface, and 
short-scale convection (with a finite wave-number) governed by surface tension gradients. A 
typical neutral stability curve for this system has two minima, corresponding to two critical 



Marangoni numbers [19|. The one minimum. Ma/, is at a zero wave-number and corresponds to 
a long-scale instability driven by deformation of the free surface coupled with surface tension 
gradients. The other minimum, Ma^, is at a non-zero wave-number, and corresponds to a 



^This type of convection can be also driven by mass rather than heat transfer |jl| 
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short-scale instability caused by surface tension gradients only (surface deformations are not 
important in this regime). 

When the two thresholds are close to each other, the nonlinear interaction between the two 
modes of Marangoni instability with different scales can lead to various secondary instabilities 
of the short-scale convection generating various stationary and temporally modulated patterns 

|2^]. The nonlinear evolution and interaction of the two modes is described by a system of 
coupled equations which can be written in the following scaled form 

At = A + A^^-\A\'^A + Ah, (2) 

ht = -mh^x - wh^xxx + sl^lL- 

Here A is the amplitude of the short-scale convection, h is the deformation of the free surface, 
X and t are long-scale space coordinate and slow time. The form of the equation for h is 
determined by the fact that h is conserved due to mass conservation. In the non-conserved case 
certain aspects of the interaction between short- and long-scale structures has been discussed 



in [21|. The parameter m is proportional to the difference between the two thresholds, Ma^ — 
Ma^, w > {) characterizes the effect of the Laplace pressure, and s > is the computable 
coupling parameter which was shown to be positive 0. Note that for the terms in the second 
equation in (^) to be of the same order, the capillary number, which parameter w is proportional 
to, has to be large (which is indeed the case in most of experimental conditions). The coupling 
between A and h reflects the following mechanism of the modes interaction: as expressed by the 
coupling term in the first equation, the short-scale convection is more intensive beneath surface 
elevations; since s > 0, the convection, in turn, suppresses the long-scale deformational mode 
by changing the average surface temperature or concentration Q. This stabilization effect of 
the short-scale convection mode allows the difference between the critical Marangoni numbers 
Ma; and Ma^ to be 0{\) and still the deformational instability to be restricted to long scales 
and have small growth rate. However, this stabilization fails if Ma^— Ma/ is too large (see [pO| ). 

In Ref.||l|, the effect of the mode interaction on the stability of the short-scale convection 
was studied only in the vicinity of the band center (yl = 1, /i = 0). In the present paper we 
extend the analysis to the full wave-number band of the short-scale convection and show that 
the interaction of the two modes strongly affects the width of the band and can even lead to a 
change of the type of the instability. 



2 Linear Stability Analysis of Short-Scale Convection Rolls 

We consider the effect of surface deformations on the stability of short-scale convection rolls 
with wave-numbers lying within an 0{e) range at small supercriticality (Ma— Mas)/Mas = e^. 
This convection pattern is described by the following stationary solution of eq.(^): 

^o = yl^e*^^ /i = 0, (3) 
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where g < 1 is the dimensionless (in units of the characteristic width e of the band of ex- 
cited wave-numbers) deviation from the critical wave-number corresponding to the threshold 
Ma=Mas of the short-scale instability. 

Consider a perturbed solution in the form 

A = Aoil + ae''^+"'^ +be''''^-'^''), (4) 

where asterisks denote complex conjugates. Substituting ^ in system (^) yields the following 
implicit dispersion relation for the side-band instability: 

+ acr^ + /3cr + 7 = 0, (5) 

where 

a = —mK'^ + WK^ + 2{l — q'^ + K^), 

ry = (_m,K2-K'u;K^)(K^ + 2^2(1 -3g2))+2sK''(l-g2)^ 

We see that the interaction of the two modes leads to a cubic equation for the eigenvalues 
defining the limits of the side-band instability (similar, say, to the case of the Eckhaus instability 



of a hexagonal pattern when there is a coupling between three waves [22, ^). In our case, we 
have the coupling between two side-band modes, which correspond to longitudinal compression- 
dilatation waves propagating in opposite directions if a is complex, and the surface deformation. 

According to the Routh-Hurwitz criterion, it follows from (^ that the stationary solution 
(^ is stable if 

a > 0, 7 > 0, a/?-7 > 0. (6) 

The analysis of conditions @ shows that, depending on the values of the parameters 
m, w and s, the side-band instability in this system can be either long-wave or short-wave, 
monotonic or oscillatory, and the range of q where the solution (|3|) is stable can be either wider 
or narrower than the standard Eckhaus interval defined by (|l]). 

A detailed analysis of various possibilities is very complicated since the stability of the 
system is governed by four parameters, m, w, s, and q. We shall dwell therefore mostly on the 
case when the solution (^) becomes unstable with respect to perturbations with infinitesimal 
wave-number. The corresponding stability limits can be determined from (^) in the limit k — > 0. 
It is, however, more instructive to obtain this result using the so-called phase approximation 
which is based on the observation that in the long-wave limit the evolution of the amplitude 



of the convection pattern is slaved to the evolution of its phase [24, 25 [. Thus, we set A = 
{po + ep) exp {iqx + (j)), h = 0(e), dx = 0(e), dt = 0{e'^), and obtain from eq.dll) for the 
amplitude correction p 

P = -—(px + Po = Jl-q^- (7) 
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In our case, the amplitude near the instabihty threshold is slaved to both the phase (/> and the 
surface deformation h. The evolution of these two variables is described by the following system 
of coupled equations: 

, 1-3^2 



ht 



(s — m)h 



+ 



1 



2 hx 1 



Introducing the wave-number deviation, k i 
coupled diffusion equations: 

at 

where 



'2(lS<f)xxx- 

X, one can rewrite system 

— S 1 



(8) 



in the form of the 



(9) 
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The system is stable if all eigenvalues of matrix S have negative real parts. The characteristic 
equation leads to the following stability conditions: 



3g^-l + (m-s)(l-g^) <0, 
{3q^ - l){m- s) + 2q'^s > 0. 



(10) 
(11) 



It should be noted that this system obtained in the long-wave limit does not give the 
complete stability information. It can be shown that there always exists an interval of m where 
the relevant instability is a short-wave instability that can be either monotonic or oscillatory. 

We shall dwell on the case when the instability is long- wave in the band center, and no 
short-wave oscillatory instability arises. A typical stability diagram is shown in Fig.l for s = 
3, w = 15. Fig.l is true provided {s < 2w, s < 1 + ^/Sw} (see Ref.[||) and w is larger than a 
certain value w*{s), respectively (for s = 3, w > w* S). Fig.l presents the region in the plane 
(m, q^) in which the solution (^) is stable. This region is bounded by the lines AB, BC, CD and 
DE. The lines AB and DE are the loci of q^ = (m — s) /{3m — s) and correspond to a monotonic 
long-wave instability (cf. ([Tl|)), while the line CD is the locus oi q'^ = {1 + s — m)/{3 + s — m) 
and corresponds to an oscillatory long- wave instability (cf. ([ToD). The line BC corresponds to 
a monotonic short-wave instability and can be found analytically from the condition 7 = (see 
Eq-dl)) (it is, however, too cumbersome to be reproduced here). The horizontal dashed line 
corresponds to the threshold of the standard Eckhaus instability, = 1/3. 

The coordinates of the points A, B, C, D, and E can be found analytically: A— > {— oo, 1/3}, 
B={(s-s^)/6, (5s+s^)/(3(s+s^))}, s^ = Vs^ + ^8sw, D= {^(1+^1 + 4/s)/2, {^1 + 4/s- 



l)/(3-\/l -|- 4/s -|- 1)}, E= {s,0}; the expressions for the coordinates of the point C are very 
lengthy (for s = 3, u; = 15 in Fig.l, C= {1.97, 0.5}). 

Thus, due to the coupling between the phase evolution and the surface deformation, both 
the interval of wave-numbers corresponding to stable stationary convection described by eq. (|3|) 
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and the type of its instability depend on the parameter m, i.e. on the difference between 
the two thresholds of the Marangoni instability, Ma^— Ma/. There are four particular values 
of m {mi{s,w), m2{s,w), m^^s), m4(s)), corresponding to the points B, C, E, D in Fig.l, 
respectively, which define four regions where qualitatively different instabilities are observed. 
These regions in the (s, m)-parameter plane are shown in Fig. 2 for w = 15. 

If m < mi{s,w), solution (^) is stable when < [m — s)/{2>m — s). Since s > in our 
case, the stability interval for q is larger than in the case of the ordinary Eckhaus instability, 
q^ < 1/3; part of the energy of phase perturbations is transferred to the damped long-scale 
deformational mode. This suppresses the Eckhaus instability and makes the stability interval 
for q larger. Outside this interval of q, solution (^) is unstable and the instability is monotonic 
and long-wave, like in the standard case with no surface deformation. In the limit m —oo, 
surface deformations are damped very strongly, and the Eckhaus interval tends to its ordinary 
value. 

When mi{s, w) < m < m2{s, w), the stability threshold for q is determined by a cumbersome 
function of w and s which we do not represent here. This threshold is still larger than 1/3. The 
instability in this case is also monotonic, but short-wave. At the point B the modulation wave 
number k goes to zero and the short-wave instability continuously changes over to a long-wave 
instability. This shows that the short-wave instability is not due to a subharmonic resonance 
as it has been found, e.g., in Taylor vortex flow [p^. 

When m2{s,w) < m < m3{s) = s, the stability interval for q is given by < (1 + s — 
m)/(3 + s — m). The instability in this case is oscillatory and long-wave and can lead to slow 
oscillations of the wave-number of the short-scale roll pattern coupled with oscillations of the 
surface deformation (see sec. 3. 2 below). 

At m = m3(s) = s, the stability region coincides with the ordinary Eckhaus interval, 
q^ < 1/3. 

If m > m3(s) = s, the roll pattern corresponding to the band center, q = 0, becomes 
unstable. This coincides with the results in [|l|. Nevertheless, as long as ^23(5) < m < 7714(5), a 
stable pattern still exists under these conditions, but with a modified wave-number lying in the 
interval (m — s) /(3m — s) < < (1 -|- s — m) /(3 -|- s — m). The lower boundary of this interval 
corresponds to the long-wave monotonic instability (i.e. the usual Eckhaus instability) while 
the upper one to the long-wave oscillatory instability. This splitting of the stable band into two 
sub-bands is reminiscent of a situation encountered in parametrically excited standing waves: 
under certain conditions the wave numbers in the center of the band become Eckhaus-unstable 
while off-center wave numbers remain stable. Such a splitting was found to lead to the formation 
of stable kinks in the wave number separating domains with different wave numbers |27|. The 
nonlinear evolution of the present system in this regime will be discussed briefly below. 

For m > 7714(5), solution (^) is unstable at any wave number due to the deformational 
instability of the surface. 
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Thus, depending on m, i.e. on the difference between the threshold Ma^ of the short-scale 
convection and that of the long-scale deformational instability (Ma^), the band of stable wave 
numbers can take on two qualitatively different forms. For m < s the band is contiguous 
(hashed area in Fig. 3a) while for s < m < 7714(5) convection is stable within two disconnected 
wave- number bands (hashed area in Fig. 3b). 

In conclusion of this section let us note that if the parameter w is not large enough, the 
oscillatory short-wave side-band instability comes into play. In this case, the stability region 
shown in Fig.l slightly changes: there appears a new boundary corresponding to this short- 
wave oscillatory side-band instability which cuts off the corner at the point C and connects 
two boundaries BC and CD corresponding to the short-wave monotonic and to the long-scale 
oscillatory side-band instabilities, respectively. 



3 Nonlinear evolution of the long-wave instabilities 

In this section we shall study the nonlinear evolution of the long-wave instability near the 
thresholds given by the curves AB, CD and DE in Fig.l. 



3.1 Eckhaus instability 

The monotonic long- wave instability threshold corresponds to the curves AB and DE in Fig.l 
and is given by 

^ (12) 



"™ 3m - s 

We use the following scaling: X = ex, T = e^i, and apply the following ansatz: 

A = pexp i(p, (13) 
p = Po + eVi {X, T) + e\X, T)p2 + . . . , 
ct> = qx + £01 (X, T) + e202(X, T) + . . . , 
h = e^hi{X,T) + €^h2{X,T) + . . . , 

q = qm + (?q2 + ■■■■ 

Substituting ( |T3|) in system (^), one obtains in the first order that the change in the local 
wave-number, ki = 5x</'i, is proportional to the surface deformation /ii. 

Proceeding to the sixth order, we obtain the following nonlinear evolution equation for the 
surface deformation hi: 

dxhi + aidxxhi + cr2dxxxxhi + cr^dxx{hl) = 0, (15) 
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where 

(4m — s){s — 3m) —3m? + ms + Asw 

m(s — 3m)(5,s — 6m) 
Ss{—m'^ + ms + s) 

The evolution equation for the long-wave modulation of the wave-number ki can be obtained 
from eq.([T5|) by a simple renormalization of the coefficient (T3 a^2sqm/{s — m) (of. (0)). 
When m — > — 00, the surface is undeformable, and the coefficients in the equation for ki tend 
to those of the nonlinear evolution equation describing an ordinary Eckhaus instability in the 
real GL equation. 

It can be shown that the coefficients (T2 and (T3 are positive (in the parameter region where 
the instability is long- wave), and the instability is subcritical. Therefore, all periodic solutions 



of (15) blow up beyond the stability limit cji = indicating a break-down of (15). The same 
happens in the case of the monotonic short-scale instability, corresponding to the curve BC in 
Fig.l; our computations show that it is also subcritical. However, numerical simulations of the 
full amplitude equations (^) show that different physical processes occur when crossing the lines 
AC and DE in Fig.l, respectively. This is related to the fact that along AC the deformational 
mode is linearly damped (m < s), whereas along DE it is excited (m > s). Figs. 4, 5 show typical 
results of numerical simulations of (^. 

Fig. 4 shows the evolution of the pattern occurring when the line AC is crossed. While 
the surface deformation and the amplitude of the short-scale convection remain regular for all 
times, the wave-number goes through a true singularity. It signifies a phase slip, which reduces 
the total phase in the system by 27r. Following this phase slip, the pattern relaxes to a periodic 
pattern with a wave-number in the stable band. 

The evolution below the line DE is shown in Fig. 5. Here the deformational mode is excited; 



it has been shown previously that the instability is subcritical at q = [20|. The simulation 
of the coupled amplitude equations (|2|) indicates that the instability does indeed not saturate 
which results in a singularity in the surface deformation and in the amplitude of the short-scale 
convection. Thus, even the amplitude equations (^) become invalid. This result suggests a 
dry-out of the fluid film as has been observed experimentally p^ ]. 

Despite the splitting of the wave-number band, no stable wave-number kinks could be 



found. This indicates that the higher-order contributions to the phase equation (15), which 
were not calculated, differ from those arising solely from a non-monotonicity of the phase 



diffusion coefficient |27|. 

It should be noted that at the point D the coefficients ai, o"2 and 173 tend to infinity, while 
at the point E the coefficient cJi vanishes regardless of the value of q2', the system dynamics 
in the vicinity of these points is more complex; it cannot be described by eq.(lT5|) and requires 
separate analysis. 
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3.2 Oscillatory long-wave instability 

Along line CD in Fig.l an oscillatory long- wave instability arises. Its threshold is given by 

2 1 + s — m 
6 + s — m 

To study its nonlinear behavior we set A = pexpicj), introduce the slow coordinates X 
Tq = e^t, T = e^t, and expand 

P = Po + mix, To,T) + e\X, To,T)p2 + e'^iX, To,T)p^ + ..., 
ct> = qx + 00 (X, To, T) + e</.i(X, Tq, T) + e^2{X, Tq, T) + e^^{X, Tq, T) + . . . , 
h = ehi{X, To, T) + e^h2{X, Tq, T) + e^^h^iX, Tq, T) + . . . , 

q = qo + (?q2 + • • • , 



(16) 



ex, 



(17) 



with 



1 



Substituting (|17|) in system (y), we obtain in the first order a linear system for hi and = dx(t>o 
which, after appropriate diagonalization, can be reduced to a linear free Schrodinger equation 



(18) 



where 



V' = ((m - s + iuoQ)hi - kQ)/{2iujQ), uoq = \J -m? + s{m + 1). 

Eq.([T^) describes waves propagating along the convection pattern and modulating its wave- 
number and the surface deformation. The phase velocity of the waves Vp depends on the wave 
number k of the modulation, as Vp = ujqk. The wave-number /cq and the surface deformation 
hi are proportional to each other 

s — m + ibjQ 



ko = 6hi, 6 



2sqo 



(19) 



and there is a phase shift between oscillations of the wave-number and the surface deformation 
equal to arg(5). 

At the onset of the oscillatory long- wave instability, a whole band of wave-numbers ~ 0(e) 



is excited. However, a strong (quadratic) dispersion of the waves described by (18) allows only 
quartic interactions between waves with different wave-numbers within the excited interval. The 
evolution of the wave amplitudes can thus be described in the third order of the perturbation 
theory by a system of Landau equations. A similar situation occured in a study of long-scale 



oscillatory double-diffusive convection ||29 |. 

In order to study the wave interaction, we consider two pairs of oppositely propagating 
waves having wave- numbers ki and K2, and present ko and hi in the following form: 



yhij 



Ai{T)e'^^ +Bi{T)e'"^ + A2{T)e'''i + B2{T)e 



je 



id:: 



J9n 



+ c.c, 



(20) 
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^1,2 — '^0^1,2^0 ± Ki^2X, 

where the ampUtudes Ai 2 and Bi 2 undergo slow evolution on the time scale T = e'^t. 



Using ( pol) and (17), we obtain from (g), as solvability conditions in the fourth and in the 
fifth orders, the following system of Landau equations for the amplitudes Ai^2 and -61,2: 

ii = 71A1 + A^il^ipAi + X\i\Bi\^Ai + A?2l^2|'^i + A?2|52p^i, (21) 
Bi = -fiBi + A^il^ip^i + X\Mi\^Bi + Xt2\B2\^Bi + aJsI^sI^^i, 

'^22l-D2| -02-1-^221^21 -D2 i- -^21 1 -^l I -D2 i- ^21 1 ^1 1 -^2- 



A2 = I2A2 + A^2l^2p^2 + A^2l^2|^^2 + \2l\Ai\'^A2 + \\i\Bi\'^ A2 
B2 = I2B2 + \22\B2f' B2 + A22|^2p-B2 + A21I-B1P-B2 -|- A2i|^lP-B2 



The linear growth rate coefficients 7^ and the Landau coefficients A^'^ (n = 1,2) have the 
following form: 

-Jn = q2a{'^^s)Kl- (3{m,s,w)Kn, = Xa,bim,s)Kl, n = l,2, (22) 

where K(d;(m, s)) > 0, ^{[3{m, s)) > if m2 < m < 771,4. The Landau coefficients A^^ (n 7^ p) 
are complicated functions of rri, s, ki, K2 satisfying the relationships 

Aia^Ki, K2) = Aaf (k2, ^i), X\2{ki, ki) = A2i(ki, ki) = A^i, 
A?2('^i. '^i) = A2i(ki, Ki) = f{m, s)A?i. 



The real parts 5R(Aa_b(m, s)) of the functions Aa^b(m,s) in Eq.(|22D are shown in Fig.6 as 



functions of m. They become unbounded in the vicinity of the point D (m = s(l + i/l + 4/s) /2) 



where system (21) is not valid. 

System (^Tj) allows to study the competition between travelling and standing waves. Let 
us consider two waves with a fixed wave-number propagating in opposite directions. Setting 



^2 = -B2 = in (21), we study the stability of the two solutions of system (|2l| ) 

1) Aj = --^,51 = ^2 = ^2 = (travelling wave); (23) 

2) Al = bI = b ^ ^2 = ^2 = (standing wave); (24) 

with respect to infinitesimal perturbations of the amplitudes Ai and Bi. The symbols of the 
Landau coefficients in (^), ( p^ and below denote now only the real parts of the respective 
coefficients. 

The stability analysis shows that 

a) travelling waves are stable if 

A?i < 0, X\, < 0, A?i/A5i < 1; (25) 
h) standing waves are stable if 

A?i<0, \xiJxU>i; (26) 
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otherwise, the system undergoes a subcritical instabihty which may lead to a blow-up and may 
return the system back to the stable region. We have performed some numerical simulations of 
(|2|) which support this expectation. 

We have found (see Fig. 6) that in our case {w = 15, ,3 = 3), when 2.135 < m < 3.361 
standing waves are stable, and when 3.361 < m < 3.745 traveling waves are stable. In the rest 
intervals of m within the interval m2 < m < rui^ close to the points C and D, both types of 
waves are subcritically unstable. 

Let us now study the stability of standing and travelling waves with respect to pertur- 
bations with different wave-numbers, i.e. the stability of solutions (|23|), (p^ ) with respect to 
perturbations of the amplitudes A2 and B2 ■ The stability analysis shows that 

a) travelling waves are stable if, besides (|2^), the following conditions are fulfilled: 

72-71^ <0, 72-71^ <0; (27) 



b) standing waves are stable if, besides (26), the following condition holds: 



72-7ia^<0; (28) 



-^21 + 
All + -^11 

The analysis of conditions ( p7| ) and (^) shows that within the intervals defined by ( p5[ ) 
and ( p6[) there exist intervals of wave-numbers where the waves are also stable to perturbations 
with an arbitrary wave-number. Fig. 7 shows how these intervals change with the variation of 
the parameter m along curve CD in Fig.l. It can be seen that standing waves are stable at 
intermediate wave-numbers, while travelling waves can be stable at infinitesimal wave-number 
as well. 

Thus, the oscillatory long- wave Eckhaus instability in Marangoni convection with de- 
formable interface can generate either standing or travelling compression-dilatation waves prop- 
agating along the periodic pattern and modulating its wavelength together with the deformation 
of the free surface. Typical results of numerical simulations of @j are given in Fig. 8a, b. It shows 
space-time diagrams for typical standing waves and traveling waves, respectively. 
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Figure 1. The dependence of the wavenumber q corresponding to the onset of the various 
side-band instabihties on the parameter m (w = 15, s = 3). Different sections of the sohd hne 
correspond to different types of the instabihty: AB - long- wave monotonic (Eckhaus); BC - 
short-wave monotonic; CD - long- wave oscillatory; DE - long- wave monotonic (Eckhaus). The 
dashed line presents the standard Eckhaus instability threshold q"^ = 1/3. 

Figure 2. Regions in the (m, s)-parameter plane corresponding to different types of the 
side-band instability (at w = 15): ML - monotonic long- wave (Eckhaus), MS - monotonic 
short-wave, OL - oscillatory long-wave. 

Figure 3. Two qualitatively different types of stable wave-number bands for the short-scale 
Marangoni convection near the instability threshold: 

a) m < m^: the stable region (hashed) is larger than in the standard case (dashed line), and 
its boundary (solid line) can correspond to either a monotonic, long- or short-wave instability, 
or to an oscillatory long-wave instability; 

b) 771-3 < rn < 1714: the band center is unstable, and the stable region (hashed) is split up into 
two parts; its inner boundary corresponds to a monotonic long-wave instability and the outer 
one to a long-wave oscillatory instability. The dotted line is the neutral stability curve. 

Figure 4. Typical space-time diagrams for the evolution of the long-wave monotonic instability 
along line AC (see Fig.l) as obtained by numerical simulation of (^): the surface deformation 
a) and the amplitude of the short-scale convection b) remain regular while the wavenumber 
passes through a singularity indicating the phase slip. 

Figure 5. Typical space-time diagrams for the evolution of the long-wave monotonic insta- 
bility along line DE (see Fig.l) as obtained by numerical simulation of (^). Here the surface 
deformation a) and the amplitude of the short-scale convection b) become singular indicating 
dry-out of the fluid film. 

Figure 6. The real parts of the interaction coefficients 3?(Aa(m)) (curve a) and 3?(Af,(m)) 
(curve b) from Eq.(p2[) at s = 3, w = 15, as functions of m. The vertical asymptote corre- 
sponds to m = 7714(5) (point D in Fig.l). 

Figure 7. Regions of stability of standing and travelling waves with a wavenumber ki (lying 
within the range of linearly excited modes) with respect to perturbations with a wavenumber 
kappa2 (lying within the same region), at different values of the parameter ttt,: s - stable, u - 
unstable. 
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Figure 8. a) Typical space-time diagram of a standing wave arising from the long-wave 
oscillatory side-band instability of the periodic short-scale convection. 

b) Typical space-time diagram of a traveling wave arising from the long-wave oscillatory side- 
band instability of the periodic short-scale convection. 
Both diagrams were obtained by numerical simulation of (^). 
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